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Abstract 

It is widely believed that quadratic divergences severely restrict natural con- 
structions of particle physics models beyond the standard model (SM). Super- 
symmetry provides a beautiful solution, but the recent LHC experiments have 
excluded large parameter regions of supersymmetric extensions of the SM. It 
will now be important to reconsider whether we have been misinterpreting the 
quadratic divergences in field theories. In this paper, we revisit the problem 
from the viewpoint of the Wilsonian renormalization group and argue that 
quadratic divergences, which can always be absorbed into a position of the 
critical surface, should be simply subtracted in model constructions. Such a 
picture gives another justification to the argument [5] that the scale invariance 
of the SM, except for the soft-breaking terms, is an alternative solution to the 
naturalness problem. It also largely broadens possibilities of model construc- 
tions beyond the SM since we just need to take care of logarithmic divergences, 
which cause mixings of various physical scales and runnings of couplings. 
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1 Introduction 



The hierarchy problem pQ , the stability of the weak scale against Planck or GUT 
scales, is considered to be an important guiding principle to construct a model 
beyond the standard model (SM). Supersymmetry is a beautiful solution, but 
the recent LHC experiments have already excluded low energy supersymmetry 
and we need to solve the little hierarchy problem as well as the \x problem. Faced 
with these difficulties, it will be important and timely to reconsider the hierarchy 
problem and ask whether we have been misinterpreting the divergences in field 
theories. 

The hierarchy problem has many faces and it is important to distinguish 
the following two types. The first one, which is most commonly referred, is why 
a scalar field can be mu ch lighter than the cutoff S ca!e0 The mass of a scalar 
field receives large radiative corrections due to quadratic divergences, and the 
tree-level mass of the Higgs field and the loop contributions of the cutoff order 
must cancel to a very high precision in order of the weak scale. Another type of 
the hierarchy problem, which is caused by logarithmic divergences, arises when 
a theory includes multiple physical scales, e.g., the weak scale and the GUT 
scale [2]. Even if the quadratic divergence is disposed of, the lower mass scale 
generically receives large radiative corrections of the higher mass scales, and a 
fine-tuning is necessary to keep the separation of the multiple scales. 

The first type of the hierarchy problem can be regarded as an academic ques- 
tion since a subtraction of the quadratic divergences is always possible without 
any physical effect on low energy dynamics o It is quite different from the 
logarithmic divergences, which play important physical roles as beta functions 
or conformal anomalies. In some formulations, one can sidestep the quadratic 
divergences and accordingly the associated hierarchy problem. A well-known 
example is to use the dimensional regularization [3]. Another way may be just 
to subtract them [JJ. A crucial feature in those formulations is that one can 
separate the subtr active and multiplicative renormalization procedures. Then 
the first type of the hierarchy problem is reduced to the naturalness of such a 
subtraction. 

As argued in ref. [5], scale invariance can be used to justify a subtracted 

lr The meaning of a cutoff is twofold. It is either a cutoff in a UV complete theory or a cutoff 
in an effective theory at which new degrees of freedom appear. We use it in both meanings 
and explain their differences depending on the situations. 

2 In a UV complete theory, quadratic divergence is only an artifact of a regularization 
procedure while, in an effective theory, it is interpreted as a boundary condition between a 
UV complete and an effective theory. In both cases, quadratic divergences can be subtracted 
without any physical effect on low energy dynamics. 
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theory, analogously to imposing the super symmetry. At the classical level, the 
SM is scale invariant except for the Higgs mass term, and a vanishing of the 
mass term would increase the symmetry of the SM. The common wisdom is that 
such increase of the symmetry cannot play any role to control the divergences 
because scale invariance is broken by logarithmic runnings of coupling constants. 
However, quadratic divergences are independent of the logarithmic divergences 
and merely an artifact of regularization procedures, and hence should be simply 
subtracted. Then the trace of the energy-momentum tensor becomes 

e£ = Am 2 H*H + p Xi O i , (1.1) 

where Am 2 is not proportional to the cutoff squared A 2 but to m 2 . The 
quadratic divergences are subtracted and the mass term gets multiplicatively 
renormalized. Here, the anomalous term as well as the mass term is regarded as 
soft-breaking terms of the scale invariance since they do not generate quadratic 
divergences. 

In this paper, we first argue that a subtraction of quadratic divergences is 
naturally performed from the Wilsonian renormalization group (RG) [6] point of 
view, and give another justification for a subtracted theory. In the Wilsonian 
RG, quadratic divergences determine a position of the critical surface in the 
theory space, and the scaling behavior of RG flows around the critical surface 
is determined only by the logarithmic divergences. The subtraction of the 
quadratic divergences can thus be performed by the position of the critical 
surface. Then, the subtraction is interpreted as a choice of the parameterization, 
i.e., a coordinate transformation, in the theory space. The fine-tuning we need 
to perform is the distance of bare parameters from the critical surface, and has 
nothing to do with the position of the critical surface itself. Hence, quadratic 
divergences are not the real issue of the fine-tuning problem [8] . 

If we consider an effective low energy theory with a finite cutoff, we en- 
counter quadratic divergences of the cutoff order. When we embed the effective 
theory in a UV complete fundamental theory, however, they are compensated 
by the same kind of divergences arising from the integrations above the cutoff 
of the effective theory. The subtraction in the effective theory is then justified 
by the boundary condition at the cutoff. Once the boundary condition is deter- 
mined by the dynamics in the UV complete fundamental theory, the quadratic 
divergences that appear in the effective theory can be legitimately subtracted 
by the boundary condition. One can also use the above arguments of coordi- 
nate transformation in the theory space to justify the subtraction within the 
effective theory. 

The second type of the hierarchy problem is more physical and should be 
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taken with much care. In the Wilsonian RG, it is formulated as a radiative mix- 
ing of multiple relevant operators, which is caused by logarithmic divergences. 
The lower mass scale is affected by higher scales through RG transformations. 
The mixing is physical and of course cannot be simply subtracted. Hence it 
gives a strong constraint on natural model building. 

The above arguments broaden possibilities of model constructions beyond 
the SM. In particular, nonsupersymmetric models with quadratic divergences 
but no large logarithmic mixings can be good candidates of models beyond the 
SM. Examples of such models are z/MSM [9] or a classically conformal TeV scale 
B-L model [TO]. 

The paper is organized as follows. We first discuss quadratic divergences in 
the Wilsonian RG in section [21 by studying a (ft 4 field theory at the one-loop 
order. We then argue that the quadratic divergence is naturally subtracted 
and is not the real issue of the fine-tuning problem in section [3j We further 
discuss logarithmic divergences and the second type of the hierarchy problem 
in section |H Finally in section [5l we show that our statements hold at all orders 
in perturbations. The last section is devoted to conclusions and discussions. 

2 RG flows of (f) A theory at one-loop 

In this section we explain the role of quadratic divergences in the Wilsonian 
renormalization group (RG). As an example, we consider a scalar field theory 
in d=4 at the one-loop approximation. The theory has both of the quadratic 
and logarithmic divergences, but the quadratic divergences can be completely 
absorbed into the position of the critical surface. In section [5] we see that it 
holds generally at all orders in perturbation expansions. 

We first consider a single scalar theory with a </> 4 interaction on a <i-dimensional 
Euclidean lattice. Its action is given by 



S 




+7| A / n^(2vr) d ^Epa)<A(pi)^ 2 )^3)^4), (2.1) 
4! y Ad ll(27r) rf ^ 



where the momentum integration is performed over the region 



A d 



{p\ - 7r < p l < IT, Vi = 1, 2, • • • 



d}. 



(2.2) 



All the quantities in the action, i.e., the parameters m 
4>, and the cutoff A = tt, are dimensionless. 



2 and A, the scalar field 
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An RG transformation can be denned by the following two steps. 
Step 1: Integration over higher momentum modes 

For simplicity, we introduce a sharp boundary at p = tt/N with a constant 
N > 1 and divide the integration region A d into two regions with lower and 
higher momenta, 

Af n = iP\~^<P l <^> Vi = l,2,. -.,<*}, (2.3) 
Alt = {P\ \P*\>^ 3i = l,2,... ,d} ■ (2.4) 

We then perform functional integrations over 4>{p) withp € A d ut . The remaining 
theory is described by an effective action for the lower momentum modes, 4>(p) 
with p G Af n . 
Step 2: Rescalings 

We then rescale the momentum p and the field 4>{p) as 

p' = Np , 

</>'(p') = N- e <P(p). (2.5) 

This rescaling of momenta makes the integration region back to the original 
one, A d . The scaling dimension 9 can be chosen so that the scalar field has the 
canonical kinetic term. (—9) is the mass dimension of the field (j)(p) and given 
by the canonical value 9 = (d + 2)/2 near the Gaussian fixed pointo 

In the d=4 (ft 4 theory, we can restrict RG transformations within the space 
of two parameters m 2 and A. The other operators, such as c/> 6 or p 4 (/) 2 , become 
less and less important when repeating the RG transformations and eventually 
become negligible: they are irrelevant operators. Hence, the restriction to con- 
sider the RG transformations within the subspace is justified. The resulting 
theory has the same form as the original one in f)2. 1 [) . but with the parameters 
changed. 

We perform the above functional integrations over 4>{p) with p £ A d ut by 
perturbative expansions with respect to the coupling A. At the one-loop order, 
the above two steps give the following changes of the parameters: 

m 2/ = N 2e ~ d {m 2 + Cl \-c 2 m 2 \) , (2.6) 
A' = N 4e ~ 3d {\ - 3c 2 A 2 ) , (2.7) 

3 If we consider RG flows near a nontrivial fixed point, 6 deviates from the canonical value 
due to an anomalous dimension of the field. Near the Gaussian fixed point, we can keep 6 
to be the canonical one and absorb all effects of the wave function renormalization in the 
coefficients c n given below. At one-loop of the <p 4 theory, a wave function renormalization is 
absent. 
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(a) 



(b) 



Figure 1: (a) Feynman diagrams that contribute to the mass renormalization 
transformation (|2.6p . The cross represents a mass insertion, (b) A diagram in 
higher order in m 2 , which does not contribute to (|2.6|) in the limit (|2.1U|) . 



where c\ and C2 are positive constants 

d d q 1 



ci 



I f d q 1 

2j A * WTw >0, ( 



out 
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The prefactors N 2e ~ d and _/V 40-M come from the rescalings (Step 2), while c n 
are contributions from the integrations of higher momentum modes (Step 1). 
The mass transformation (|2.6jl is given by the diagrams in Figure Ufa) • Since 
the integration is performed in the UV region A d ut and no IR divergences occur, 
we can expand the propagator with respect to m 2 . In eq. (|2.6|) we took the first 
two terms in the expansion. Higher order terms in m 2 , such as c^m^X which 
comes from a diagram in Figure QJ b), are highly suppressed since we suppose 
that 

2 

m 

JT « 1 • (2-10) 

Similarly the coupling transformation (12. 7p at one-loop is determined by the 
diagram in Figure [2f a). Again, higher order terms in m 2 , such as c' 2 m 2 X 2 for 
Figure Efb), do not contribute to (|2.7p due to the assumption (|2.10p . There are 
no wave function renormalizations at the one-loop order. The dependence of 
the constants c n on the cutoff A is easily derived □ By evaluating the integrals 
of (ELS) and (E2D, we find 



ci oc A d ~ 2 (l -AM d - 2 )) A 2 (l-iV- 2 ) ford = 4, (2.11) 

c 2 oc A d ~ A (l -iV-( d - 4 )) A°lniV for d = 4 . (2.12) 



4 In taking the continuum limit, we consider a limit where the physical momentum becomes 
infinitesimal compared to the cutoff A = ir. Hence, A = n corresponds to a large momentum 
scale in the unit of the physical scale. 
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(a) (b) 

Figure 2: (a) A Feynman diagram that contributes to the coupling renormal- 
ization transformation (|2.7p . (b) A diagram in higher order in ?re 2 , which does 
not contribute to ([22]) in the limit (|2TTU|) . 



Hence, c\ and c 2 reflect the quadratic and the logarithmic divergences, respec- 
tively. 

By performing the RG transformations (|2.6p and (|2,7p several times, one 
obtains RG flows in the theory space, the space spanned by the parameters A 
and m 2 in this case. Since eq. (12. 7j) depends only on A, but not on m 2 , the flow 
of the coupling constant A is determined only by using (|2.7p . For d ^ 4, the 
equation (|2.7p is rewritten as 



I - 1 = N-W-™) (\-^]+ O(X) , A* = — M 1 . (2.13) 
A' A* \X X* J w ' 3c 2 v ' 

After performing the transformation n times, one obtains 

— - — = AMMQn _ / 2 .14) 

An A* VAo A*y 

Here, A n is the renormalized coupling constant after re-times RG transforma- 
tions. There are two fixed points of the RG transformation: one with A = 0, 
and the other with A = A*. For d = 4, two fixed points coincide at A = and 
one finds 

T^ = T*- + 3c 2 n. (2.15) 

An Ao 

If Ao > 0, A„ approaches the fixed point A = as one increases n. 

In order to obtain the flow in the direction of m 2 , we rewrite eq. (|2.6p . by 
using (|2.7p . as 



m 2 ' - m 2 (A') = N 2e ~ d {l - c 2 A)(m 2 - m 2 (A)) , (2.16) 
with a function m 2 (A), determined up to this order as 

m 2 (A) = C \ la A . (2.17) 
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Performing the RG transformation n times, one obtains 

n-1 

ml-ml{\ n ) = NW-^Hil-czK) (m 2 - m 2 c (X Q )) 

i=0 

n-1 

~ iV( 2 ^"exp(-c 2 ^A t ) (mg-m^(Ao)). (2.18) 

i=0 

The equation m 2 = m 2 (X) determines the position of the critical line, and 
eq. (|2.18p shows how the distance from the critical line scales under the RG 
transformations . 

The RG flow is given by the two equations (pTTBj) and (gUj) , or (j2T5]) for 
d = 4. A schematic picture of the flow for the d = 4 case is drawn in Figure [3j 
The fixed point of the RG flow is given by m? = A = 0, i.e., the Gaussian fixed 
point. If the initial parameters are exactly on the critical line m 2 = m 2 (X), they 
continue to be there after RG transformations and approach the fixed point. 
We consider only the region A > for the stability of vacuum. A theory on 
the critical line is a massless theory. If the initial parameters are off the critical 
line, they depart from it with A decreasing under the RG transformations. 

As we stressed in the introduction, the constant c\, which reflects the 
quadratic divergences, is completely absorbed into the position of the criti- 
cal line, i.e., the definition of the function m^A) in (|2.17|) . and the scaling 
behavior (|2,18p of the RG flow around the critical line is determined only by C2, 
which reflects the logarithmic divergences. It is a very important message we 
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can read from the Wilsonian treatment of the renormalization group. As we see 
in the next section, it corresponds to the fact that in ordinary perturbative cal- 
culations, quadratic divergences can always be subtracted without any physical 
effects on the dynamics of field theories, unlike logarithmic divergences. 

3 Quadratic divergences, subtractions, and the fine- 
tuning problem 

We now interpret the subtraction and the fine-tuning problem in the framework 
of the Wilsonian RG, and argue that the quadratic divergence can be naturally 
subtracted and is not the real issue of the hierarchy problem. 

3.1 Continuum limit 

We first review how to take the continuum limit. In the constructive formulation 
of a field theory, one usually takes the continuum limit by simultaneously letting 
the parameters close to the critical surface and taking the cutoff to infinity. For 
a scalar theory in d = 4, however, one cannot construct an interacting theory 
in this way, which is well-known as the triviality or the Landau singularity 
problem. We thus consider a theory at a large but finite cutoff with bare 
parameters very close to the critical surface. 
We first write the RG equation (|2.18p as 

m 2 - m c 2 (A ) = AT-^-^E^ (m 2 n - m 2 c (X n )) . (3.1) 

The parameters ml and Ao on the left-hand side (LHS) describe the bare param- 
eters, while m\ and X n on the right-hand side (RHS) describe the renormalized 
ones. Note that all the quantities in this relation are dimensionless, including 
the cutoff A = 7r. We now introduce dimensionful cutoffs. First we introduce 
the low energy scale A n = M, e.g., M = 100 GeV. The renormalized param- 
eters m n and A n are defined at this scale. The higher scale where the bare 
parameters too, Ao are defined is given by 

A = N n A n = N n M . (3.2) 

Let us express everything in terms of dimensionful physical quantities in- 
stead of the dimensionless lattice parameters. Physical momentum p is defined 
by 

Pk = ^P, A k = N n ~ k M, (3.3) 
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for each cutoff theory at k = 0, l,--- , n. We can similarly define the field 
as 4>k(pk) = (Afc / 'A)~ ""a - ' (f){jp) . The dimensionless parameters m 2 ,A are also 
replaced by a dimensionful mass and a coupling 

In terms of these physical quantities, the RG equation (|3.1j) can be rewritten 

as 

mg - m 2 (A ) = N^- 2e+d >e c ^ ^ (™ 2 _ ^2 (An)) (3 5) 

The critical value m 2 (A) is also rescaled as 

m2(A*) = f ^ m c 2 (A fc ) , (3.6) 

which reflects quadratic divergences of the position of the critical line, as shown 
in the previous section. For the Gaussian fixed point, we have 9 = (d + 2)/2, 
and hence Ar( 2 - 2<? + d )« = 1. Eq. ([33]) becomes 

(m 2 -m 2 (A )) = e c ^r= 1 A i _^2 (An)) (3J) 

The prefactor e C2 Si=o Ai represents the logarithmic running of the mass pa- 
rameter. Indeed, from (|2.15p . one can find Aj ~ l/3c2«, and thus the prefactor 
behaves as e C2 ^=o Ai ~ ra 1 / 3 ~ (InAo) 1 / 3 . Its n-dependence is weaker than a 
power-law running behavior N xn ~ (Ao) x with some constant x. 

In ordinary perturbative calculations of field theories, we first choose a reg- 
ularization method and a renormalization prescription to deal with divergences. 
Quadratic divergences are simply subtracted with appropriate renormalization 
conditions. The so called fine-tuning problem is that we need to fine-tune the 
bare mass m 2 , of a scalar field at the cutoff scale Ao, so that the renormalized 
mass m 2 R , 

rhR = m o + hi(A /A„)) A 2 , + /3(A, m(A /A n )) m 2 , 

= (1 + /3(A, ln(A /A n ))) (m 2 + a'(A, ln(A /A re )) A 2 ) , (3.8) 

is held at the desired value. Here a and j3 are functions of the coupling A and 
ln(Ao/A n ), and assumed to be independent of the mass parameter ttiq. Namely, 
the mass-independent renormalization scheme [114112] is chosen. When we apply 
(|3.8p to the Higgs mass, this looks very unnatural because, if Ao S> rfiR = m\y 5 
the tree-level mass m 2 , of the Higgs particle and the loop contributions oAq 
must cancel to a very high precision in order of the weak scale m 2 ^. There 
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is also a subtlety in the mass-independent renormalization for theories with 
quadratic divergences [121 H] • 

Let us make a comparison between eq. (|3,7|) and eq. (|3.8|) . The bare mass 
parameter fh\ and the quadratic divergence a Kq in the RHS of (|3.8|> corre- 
spond to the bare mass and the critical value mj?(Ao) in the LHS of (|3.7p . 
respectively. The renormalized mass mfj of (|3.8p is given by — rh^(X n ) in 
the RHS of (13. 7p . The multiplicative renormalization factor (1 + /3) in (13. 8p 
corresponds to the factor for the logarithmic running e~ C2 ^=o Xi in (|3.7p . 

3.2 Subtractive renormalization in Wilsonian RG 

In the Wilsonian RG, the quadratic divergence c\ is absorbed into the position 
of the critical line rh^(X). Note also that all observable quantities like the cor- 
relation length are determined by the distance of the mass parameter rh\ from 
the critical line m^(Ao), not by the value jbq itself. The fact that the difference 
ttiq — m^(Ao) is the physically relevant quantity gives a natural interpretation 
for the subtractive renormalization of the quadratic divergences. 

As we stressed in the introduction, an important and necessary feature for 
the subtracted theories, as in [31 dj [5] , where one can sidestep the quadratic 
divergences, is that one can separate the subtractive and the multiplicative 
renormalization procedures. Looking at eq. (|3.7p . the prefactor e C2 ^<=o Ai corre- 
sponds to the multiplicative renormalization. Hence, eq. (|3.7p clearly separates 
the subtractive, TOq — m;?(Ao), and the multiplicative, e C2 ^--»=o x ', renormaliza- 
tion procedures. 

Another remarkable result from eq. (|3,7p is as follows. Neither the prefac- 
tor e C2 ^i=o Xi nor the critical value rn^(X) depends on the mass parameter mo- 
Hence, eq. (|3.7p gives an explicit realization of the mass-independent renormal- 
ization scheme [T2] . As discussed in |12^ U] , whether or not one can formulate 
a mass-independent renormalization scheme in theories with quadratic diver- 
gences was a nontrivial and subtle problem. 

Note that the subtraction is not arbitrary in the Wilsonian RG. Once we 
choose a scheme to calculate the RG flow, the critical line is given unambigu- 
ously. Seeming ambiguities only come from our lack of calculability of the exact 
RG flows, but RG flows are in principle exactly determined once we choose a 
scheme. The only fine-tuning we need is the distance from the critical surface 
to the bare parameters. 

Before discussing the tuning of the distance, let us convince ourselves that 
the position of the critical surface has nothing to do with the fine-tuning prob- 
lem by comparing two theories: one with a nonzero c±, and the other with a 
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vanishing c\ . The latter theory corresponds to a theory with vanishing quadratic 
divergences. The critical lines are given by m 2 7^ for c\ 7^ 0, and m 2 = for 
ci = 0. Is the first theory more unnatural than the second one? The feeling of 
unnaturalness simply comes from seeming ambiguities of calculating the posi- 
tion of the critical line, but as stressed, once a scheme is given it is determined 
unambiguously and there are no distinctions between these two theories. In 
other words, we can take new coordinates of the theory space such that the 
critical line is parallel to an axis of the coordinates. Then a theory with c\ 7^ 
looks the same as a theory with c\ = 0. The change of coordinates can be 
determined exactly once we choose a scheme. In the dimensional regulariza- 
tion, quadratic divergences do not appear. We may say that this corresponds 
to taking a new coordinate by m^ ew = m 2 — m 2 (X). Then the critical line is 
given by m^ cw C (A) = even for a theory with c\ 7^ 0. Therefore the subtraction 
of the quadratic divergences is simply a choice of the coordinates of the theory 
space. In this sense, the quadratic divergences are naturally subtracted and are 
not the real issue of the hierarchy problem. 

Finally, we consider the tuning of the distance from the critical line, namely, 
the fine-tuning to choose the bare mass near the critical value with the precision 
of the weak scale. Such a fine-tuning is kinematical in the following sense. 
Looking at eq. (|3. 1 j) . we see that this kind of tuning comes from the factor 
_/y-(20-d)n = i\r _2n for the Gaussian fixed point, but it reflects nothing but 
the canonical dimension of the mass and has nothing to do with the quadratic 
divergences. This fine-tuning is necessary as long as the dimension of mass 
square is two. The prefactor e C2 ^ Xi simply gives a logarithmic scaling factor 
ra 1 / 3 , and does not change the situation much. 

If we consider a nontrivial fixed point instead of the Gaussian one, the 
canonical scaling dimension is modified by anomalous dimensions of the mass 
and the field. For example, if 8 = {d + 2)/2 — 8 and the coupling constant at 
the fixed point is given by A*, the relation (|3,7p is modified aaj 

{ml - m 2 c (X )) ~ {N 25 e c > x *T(m 2 n - m c 2 (A n )) . (3.9) 

Thus, large anomalous dimensions can relax the condition of how precisely we 
need to define the bare mass near the critical line. If one feels uncomfortable 
with the relation (|3.7p . a possible resolution will be to construct a theory with 
large anomalous dimensions. 

5 The prefactor e C2 ^=a Xi in 1|3.7[) gives not only the factor e C2A * n in (|3.9|l . but also an 
extra factor e C2 ^i=o ( A i _A *). Since it has a weaker n-dependence than a power-law behavior 
N xn ~ A 21 , we neglected it here. The wave function renormalization also gives a deviation 
from 6, 8 and an extra factor neglected in (|3.9[) . 
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4 Mixing of multiple scales 



In the previous two sections, we considered a theory with a single physical scale 
besides the cutoff scale. In this section, we study an RG flow in a theory with 
hierarchically separated multiple scales. These scales are mixed by radiative 
corrections associated with the logarithmic divergences. This causes the second 
type of the hierarchy problem. We emphasize that it again has nothing to do 
with the quadratic divergences. 

For simplicity, we consider a theory with multiple scalar fields on a d- 
dimensional Euclidean lattice, whose action is given by 



s 



t 



i 



i 



(p + m a )4> a + -X aa 

a=l 



(4.1) 



with \ a/ 3 = \p a . The index a = 1,...,S labels species of the scalar fields. 
There are S mass parameters rr? a and S(S + l)/2 coupling parameters A Q « in 
this case. The momentum integration region A d is taken as in eq. (12.2j) . 

To obtain an RG transformation, we follow the two steps in section [2j At 
the one-loop order, it becomes 



,2 / 



A' 



a: 



N 



29-d 



m 2 a + ci A c 



C2 



N 



3c 2 ^(A a/3 ) 2 



N 



48 -3d 



A a /3 — C2 ^ Aq, 7 A/3 7 — 4:C2(X a pY 



(4.2) 



(4.3) 



(4.4) 



with c\ and c 2 given in (I2.8P and (I2.9p . 

By performing the RG transformations several times, one obtains RG flows 
in the theory space, S(S + 3)/2-dimensional space spanned by the parameters 
A a /3 and rr? a . We use A to represent the set of couplings \ a p in the following. 
Since eqs. (14.31) and (I4.4p depend only on \ a p, and not on m^, the flow in \ a p 
is determined by (14. 3p and (|4.4j) . For the flow in m^,, we can rewrite eq. (14. 2p 



as 



,2 / 



e (V) = N™~ d J2(S a p ~ <*\af,){m} ~ m%{\)) , (4.5) 



where m^ a (X) are functions of \ a p and are determined by the equations 

m c 2 a (A') - N 2d ~ d ^ap ~ c 2 \ a/3 )m%(\) = N 2e - d Cl ^ A a/3 , (4.6) 
P P 
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together with (|4.3|) and (|4.4p . The solutions are given as a power series of A as 

ci 

1 _ AT2(e-rf) 



"4(A) = -— £ + ^(A 2 ) • (4.7) 



The quadratic divergence in d=4, namely ci, is again absorbed into the position 
of the critical surface = m^ a (X). 

Now let us define a symmetric mixing matrix 

(Afy.)) a/ 8 = <5 a/3 - c 2 A a ^ (fc ) ~ exp(<5 a/3 - c 2 A a/3 ( fc )) , (4.8) 

and their product 

M = M in _ 1} M {n _ 2 y ■ M {0) . (4.9) 

Here \ a p(k) represents the coupling constant \ a p after fc-times RG transforma- 
tion. By performing the transformation (|4.5|) n times, one obtains 



ml [n) - m 2 ca (X (n) ) = N^ n Y^M^ (m 2 m - m^(A (0) )) • (4.10) 

In taking the cutoff very large, one needs to fine-tune the S relevant opera- 
tors, m^oj. We rewrite the RG equation (|4.10p as 

<o) -m 2 ca (\o)) = N-W- d ^Y,( M ~%eH(n) -m%{\ n) )) , (4.11) 



where 

n-l 

(M" 1 )^ = ^ + c 2 K m + ©(A 2 ) • (4.12) 

fc=0 

The equation (|4.1ip can be written in terms of the dimensionful parameters, as 
in (13.71). as 



C(A (0) ) = £ ( M ^ 1 ) a /3 ("*?(n) " ^ 2 /3(A(„))) • (4.13) 



~ 2 ^,2 

/3 



Here we used the canonical value of # = (d + 2) /2 for the Gaussian fixed point. 

We now simplify the discussion by considering a theory with two separate 
renormalized scales, e.g., 

^l(n) - ™ci( A (n)) = m w , 

™2(n) -™e2( A (n)) = ™GUT- ( 4 -14) 

It then follows from (|4.13p that we need to fine-tune the bare mass parameters 
such that the difference of {m\ — ™ci(^o)) an d radiative corrections from the 
higher physical scale ^gut are canceled to give the weak scale: 

1 ,~i -2n « (M- 1 )^ 



™gut • (4-15) 
li 
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Unlike the subtraction of the critical mass parameter m^(A), the GUT scale 
■n^GUT is a physically observable scale and cannot be subtracted. Or, in other 
words, the second term proportional to rn 2 ^^ cann °t be absorbed into the 
position of the critical surface. Hence, unless the off-diagonal element of the 
matrix (M )i2 is suppressed, we need a fine-tuning of the bare mass (m 2 — 
m^Ao)) against the GUT scale with a high precision in order of the weak scale. 
As we can see from (|4. 12[) . in order to solve the problem, we need to suppress 
either the coefficient C2, the mutual couplings X a p, or the higher scale w-qtjt- 



5 Higher orders of perturbations 

Up to now, our discussions are based on the one-loop order calculations, but our 
statements hold at all orders of perturbations in the coupling A. In this section, 
we extend the statement in section [2] to all orders. It is also straightforward to 
extend the results in the other sections in the same way. 

In the following, by using a renormalized perturbation theory, we will show 
iteratively that the statement in section 2 does hold at all orders of perturba- 
tions in the coupling constant A. We first replace the mass parameter m 2 by 
m 2 — m 2 + m 2 in the action (12.11) . 



S 



ipY + \(m 2 - m 2 (A))</> 2 + \ m 2 c m 2 + lA0 4 



(5.1) 



where the A dependence of the position of the critical surface m 2 (X) will be 
determined later in a self-consistent way. We then perform the functional inte- 
grations over the higher momentum modes (Step 1 of the RG transformation) 
by perturbative expansions with respect to m 2 — m 2 (A) and m 2 (A), as well as A. 
Such mass expansions are legitimate since the integrations are performed only 
in the UV region and free from IR divergences. An insertion of the renormalized 
mass parameter m 2 — m 2 (X) suppresses loop integrations by a factor 1/A 2 , and 
thus always appears in a combination of 

^-rtW-of') . (5.2) 



A 2 VA 2 , 

Here, we assumed that m 2 — rn 2 (X) is of order A . It was explicitly shown at 
the one-loop order in sections [21 and will be justified later in this section in an 
iterative way at higher orders in A. We can thus neglect higher order terms in 
the expansion of m 2 — m 2 (X). On the other hand, an insertion of the critical 
mass is not suppressed by 1/A 2 since the value m 2 (X) itself is of order A 2 : 

1^ = (A). (5.3) 
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Figure 4: The Feynman diagrams that contribute to /A 2 up to order A 2 . The 
dot represents an insertion of m 2 . The first two diagrams give order A 1 contri- 
butions, while the last four give order A 2 . Note that the second diagram gives 
both order A 1 and A 2 contributions. 



Though it is suppressed by the coupling constant A, it is not suppressed by the 
cutoff A. It was already shown in (12,170 at one-loop order, and can be justified 
later in this section in an iterative way. Hence, on the contrary to (|5.2p . we 
need to take higher orders of m 2 (A) in the perturbative calculations. 

As a result of the above arguments, we can schematically write the effective 
action after the integrations as 



A/N 



/(A,^f) A 2 + 9 (A,^|) (m 2 -m*(A)) 



1 , m rs 

+1 ( ' ~$> 



(5.4) 



where the functions e, /, g and h are power series with respect to A and m 2 /A 2 , 
with the coefficients dependent on N. /A 2 is given by the 2-point, 1PI diagrams 
with the external momentum fixed to be zero. They are depicted in Figure H] 
up to order A 2 , where the dot represents an insertion of m 2 . g is given by 
those diagrams with a single insertion of m? — m 2 (A). They are depicted in 
Figure [5] up to order A 2 , where the cross represents an insertion of m 2 — m 2 . 
e — 1 is given by the diagrams of /, but the first term in the expansion of 
the external momentum. They begin with the third diagram in Figure 0] at 
order A 2 . Similarly, h is given by the 4-point, 1PI diagrams with the external 
momentum fixed to be zero. They are depicted in Figure [6] up to order A 3 , 

We next perform the wave function renormalization, 4> — > e -1//2 <^>, and the 
rescaling of momenta and fields. It is Step 2 of the RG transformation, and we 
have 



1 



+ 



1 



N 2B-d e - 



1 



[f A 2 + <7(m 2 -m 2 (A))] + -N 



W-3d e -2 h 
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Figure 5: The Feynman diagrams that contribute to g up to order A 2 . The 
cross and the dot represent an insertion of m 2 — m 2 and m 2 , respectively. The 
first diagram gives an order A contribution, the second A , and the last four 
A 2 . 




Figure 6: The Feynman diagrams that contribute to h up to order A 3 . The dot 
represents an insertion of m 2 . The first diagram gives an order A 1 contribution, 
the second A 2 , and the last four A 3 . 
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We then find a generalized formula of (|2.6p and (|2.7p , as given by the equations 



N 2e-d e -i ^ A 2 + g (m 2 _ m 2(A))] , (5.5) 



A' = iV 49 - 3d e (A, ^ )- 2 fc(A, ^) • (5.6) 
As in eq. (|2.16p . we rewrite f)5.5j) as 

m 2 ' - m 2 (A') = 7V 2 ^e(A, ^)~ 1 5(A, ^) (m 2 - m 2 (A)) , (5.7) 

with 

ml(X') = iV 2 ^e(A, ^)~ 1 /(A, ^) A 2 . (5.8) 

The functional form of m 2 (A) is determined iteratively in A by solving the 
equation (|5.8p . with (|5.6p inserted in the LHS of (|5.8p . At the order of A 1 , / is 
given by the first two diagrams in Figure HI Then, the solution m 2 (A) indeed 
becomes (I2.17p . At the order of A 2 , / is given by the last four diagrams in 
Figure [H The dot in the loop of the last diagram means the order A 1 term in 
m 2 (A), while the dot in the second diagram means the order A 2 term in m 2 (A). 
e has the order A 2 term, but it does not affect (|5.8p at order A 2 . The LHS of 
(|5.8|) has both contributions from the order A' 1 term in m 2 (A') with the A 2 term 
of (15. 6p . and from the order A' 2 term in m 2 (A'). In this way, we could obtain 
the A 2 term in m 2 (A). We can also obtain higher order terms iteratively in the 
same way. 

The most important property of eq. (15. 8p is that the solution m 2 (A) is 
proportional to A 2 at all orders in perturbations. Following the same arguments 
given in section [21 we can see that m 2 (A) gives the position of the critical line, 
and (|5.6p and (|5.7p determine the scaling behavior of the RG flows around the 
critical line. We therefore find that at all orders in perturbative expansions in A, 
the quadratic divergences are completely absorbed in the position of the critical 
line, and they do not play any role in the dynamics of field theories. The RG 
flow around the critical line is determined only by the logarithmic divergences. 
Hence, the assumptions (|5.2|) and (|5.3p we adopted in this section are justified 
iteratively in A. 

As we mentioned in section [3j eq. (|5.7p means that we can separate the 
subtractive and multiplicative renormalization procedures. Eqs. (|5.6p and (|5.7p 
also give an explicit representation of the mass-independent renormalization 
scheme. The separability of quadratic and logarithmic divergences relies on 
several properties of loop integrations in Wilsonian RG. First of all, they are free 
from IR divergences and hence, we can perform the mass expansions. Secondly, 
the loop integrations are performed from A/N to A, and thus a divergence 
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from the subdiagram is compensated by the remaining part of the diagram 
with negative power of A. For instance, in the 5th diagram of Figure while 
the upper subdiagram gives a A 2 contribution, the lower one gives A -2 , which 
makes the overall diagram of order A and insensitive to the divergence from the 
subdiagram. We thus need to take care of only the overall superficial degrees 
of divergencej§. That is why the naive dimensional analysis is possible. 

6 Conclusions and discussions 

In this paper, we revisited the hierarchy problem, i.e., stability of mass of a 
scalar field against large radiative corrections, from the Wilsonian RG point of 
view. We first saw that quadratic divergences can be absorbed into a position 
of the critical surface m 2 (A), and the scaling behavior of RG flows around the 
critical surface is determined only by logarithmic divergences. The subtraction 
of the quadratic divergences is unambiguously fixed by the critical surface. In 
another word, the subtraction is interpreted as taking a new coordinate of the 
space of parameters such that m 2 ew = m 2 — m 2 (A). These arguments gave 
a natural interpretation for the subtractions, and another justification for the 
subtracted theories as in [31 [H [5]. The fine-tuning problem, i.e., the hierarchy 
between the physical scalar mass and the cutoff scale, is then reduced to a 
problem of taking the bare mass parameter close to the critical surface in taking 
the continuum limit. It has nothing to do with the quadratic divergences in 
the theory. Therefore the quadratic divergences are not the real issue of the 
hierarchy problem. If we are considering a low energy effective theory with 
an effective cutoff, the subtraction of the quadratic divergences corresponds to 
taking a boundary condition at the effective cutoff scale. Hence it has nothing 
to do with the dynamics at a lower energy scale, and when such divergences 
appear in radiative corrections, we can simply subtract them. 

We also considered another type of the hierarchy problem. If a theory con- 
sists of multiple physical scales, e.g., the weak scale m\y and the GUT scale 
m GUT besides the cutoff scale A, a mass of the lower scale mw receives large 
radiative corrections dm 2 ^ oc mQ UT log A/ 167T 2 through the logarithmic diver- 
gences. Such a mixing of physical mass scales is interpreted as a mixing of 

6 Our renormalized perturbation is reminiscent of the usual one of BPHZ. In our case, 
however, the cancellation between the loop contribution and the counter term is not exact. 
For example, at order A , the first and the second diagram in Figure[4]did not cancel in (|5.8[1 . 
It follows that, at order A 2 , the loop contribution of the upper subdiagram in the 4th diagram 
is not canceled by the dot in the 5th diagram. One can also see the same situation in the last 
two diagrams in Figure [5] However, divergences from subdiagrams cause no problem in our 
case. 
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relevant operators along the RG flows. Unlike the first type of the hierar- 
chy problem, the mass of the larger scale togut cannot simply be disposed 
of by a subtraction. In order to solve such a mixing problem, we need to 
suppress the mixing by some additional conditions. Of course, if these two 
scalars are extremely weakly coupled, the mixing can be suppressed. If the 
couplings are not so weak, we need to cancel the mixing by symmetries or 
some nontrivial dynamics. A well-known example is the supersymmetry, where 
the non-renormalization theorem assures the absence of such mixings unless 
supersymmetry is broken. 

Let us comment on scheme dependence of the subtraction. In this paper, 
we fixed one scheme to perform RG transformations. Then the critical surface 
is unambiguously determined. If we change the scheme, e.g., from the sharp cut 
off (j2.4|) of the higher mode integrations to another one [13] , the RG transforma- 
tions are changed, and so accordingly is the position of the critical surface. But 
the definition of bare parameters is correlated with the choice of the scheme. 
Hence a shift of the position of the critical surface by a change of scheme does 
not mean an ambiguity of the critical surface. Rather it corresponds to chang- 
ing coordinates of the theory space. In this sense, the subtraction of a position 
of the critical surface from the bare mass is performed for each fixed scheme 
without any ambiguity. 

What is the meaning of the coefficients of various terms in the bare action? 
In the investigation of the RG flows, we encountered two kinds of quantities, 
the mass parameter m 2 and the subtracted mass (m 2 - m c (A) 2 ). The issue of 
the fine-tuning problem, i.e., the stability against the quadratic divergences, is 
related to which quantity we should consider to be a physical parameter. In 
the renormalization-group-improved field theory, as we studied in this paper, 
the subtracted one is considered to be physical. The mass parameter itself 
depends on a choice of coordinates of the theory space, and changes scheme by 
scheme. If we want to derive low energy field theories from a more fundamental 
theory, like a string theory, the bare parameter itself, m 2 , is widely believed 
to be related to the fundamental quantities at the string scale. But since such 
a quantity is coordinate (of the theory space) dependent, we should use the 
subtracted one when we relate low energy field theories with more fundamental 
theories. The amount of subtraction is given by the boundary condition at the 
cutoff, and determined by the dynamics at higher scales in the fundamental 
theory. It is independent of the low energy dynamics. 

We are thus left with the second type of the hierarchy problem, namely a 
mixing of the weak scale with another physical scale like ttigut- We classify 
possible ways out of it: 
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1. SM up to A 



2. New physics around TeV, but nothing beyond up to A 

3. New physics at a higher scale, but extremely weakly coupled with SM 

4. New physics at a higher scale with nontrivial dynamics or symmetries 

The first possibility is to consider a model without any further physical scale 
up to the cutoff scale A. The Planck scale may play a role of a cutoff scale 
for the SM. As we saw in this paper, the quadratic divergence of the cutoff 
order can be simply subtracted and it does not cause any physical effect. In 
the second possibility, we introduce a new scale which may be coupled with the 
SM, but suppose that the new scale is not so large compared with the weak 
scale. Then even if the mixing is not so small, the weak scale does not receive 
large radiative corrections. Various kinds of TeV scale models are classified 
into this category. Some examples are z/MSM [S] and the classically conformall] 
TeV-scale B — L extended model [10] . The third one is to consider a very large 
physical scale, but with the mutual coupling suppressed to be very small. The 
final possibilities include a supersymmetric GUT, but the low energy theory of 
broken supersymmetry must be supplemented with the second type of scenario. 
If we worry about quadratic divergences, the first three categories need fine- 
tunings against the cutoff scale and are excluded by the naturalness condition. 
Hence, most model building beyond the SM has been restricted to the last 
category. Once we admit that quadratic divergence is not the real issue of the 
hierarchy problem, it broadens our possibilities of model constructions. 

The discussions of the quadratic divergences in this paper can also be ex- 
tended to the quartic divergences, and if gravity is described in terms of a 
renormalized field theory, the Wilsonian RG treatment might give a new per- 
spective of the cosmological constant problem. We hope to come back to this 
issue in future. 
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